Tunable Wigner States with Dipolar Atoms and Molecules 
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We study the few-body physics of trapped atoms or molecules with electric or magnetic dipole 
moments aligned by an external field. Using exact numerical diagonalization appropriate for the 
strongly correlated regime, as well as a classical analysis, we show how Wigner localization emerges 
with increasing coupling strength. The Wigner states exhibit nontrivial geometries due to the 
anisotropy of the interaction. This leads to transitions between different Wigner states as the tilt 
angle of the dipoles with the confining plane is changed. Intriguingly, while the individual Wigner 
states are well described by a classical analysis, the transitions between different Wigner states 
are strongly affected by quantum statistics. This can be understood by considering the interplay 
between quantum-mechanical and spatial symmetry properties. Finally, we demonstrate that our 
results are relevant to experimentally realistic systems. 



Recent experimental advances in cold quantum gases 
have placed focus on atoms or molecules with permanent 
dipole moments. A Bose-Einstein condensate of 52 Cr 
atoms has been realized [TJ |2], and recently Dy atoms 
were cooled and trapped [3]. These atom species have 
magnetic dipole moments of several Bohr magnetons. A 
promising development is the trapping and cooling of di- 
atomic molecules with electric dipole moments [4-6 J . The 
realization of a molecular fermionic 40 K 87 Rb gas was a 
significant breakthrough |7J. Dipolar gases offer access 
to a broad range of novel few- and many-body physics, 
in single traps as well as in optical lattices which has 
spurred intensive theoretical interest (see the recent re- 
views by Baranov [8 J and Lahaye et al. |9J). The attrac- 
tive part of the dipolar interaction leads to a collapse in- 
stability in three dimensions |9J . A remedy is to use traps 
of reduced dimensionality. For instance, the interaction 
between dipoles in a 2D plane is predominantly repulsive 
when they form a sufficiently large angle with respect to 
the plane (see Fig. [I]). This stabilizes the system against 
collapse [TO] . A variety of interesting many-body states 
for dipoles in 2D has been examined theoretically [T0HT2] . 

Here, we examine a 2D system of this kind in the 
regime of strong repulsive interactions using exact diago- 
nalization as well as a classical analysis. Contrary to the 
analogous Wigner states of electrons in metals [14] 

and quantum dots [15J the anisotropic dipolar interac- 
tion is shown to give rise to Wigner states with nontriv- 
ial geometries that depend on the alignment angle. This 
leads to transitions between different geometries that are 
crucially influenced by quantum statistics. We finally ar- 
gue that our results are experimentally observable. 

We consider particles with mass m and a magnetic 
(electric) dipole moment \i (d) which is aligned by an 
external field such that it lies in the xz plane, forming 
an angle © with the x axis (see Fig. [I]). The interaction 
between two electric dipoles separated by a vector r is 



V(r) 



d 2 l-3cos 2 6>. 



rd 



47re n 



(1) 



where 6 r d is the angle between a dipole moment and r. 
The particles are confined in the xy-pl&ne by a quasi-two- 




Figure 1: Dipolar particles in a quasi-2D trap in the xy plane. 
The gray circles indicate the contours of the effective interac- 
tion potentials in the plane. In the left-hand panel the dipoles 
are perpendicular (G = 90°) to the plane of motion and the 
interaction is isotropic in the plane, while in the right-hand 
panel the dipoles are tilted (G = 55°) and the interaction is 
anisotropic. 



dimensional harmonic trap Vt r &p(x,y, z) = m[ujQ(x 2 + 
y 2 ) + u 2 z 2 ]/2, with u z ^> uq so that the particles are 
in the ground state orbital of the z direction. The cor- 
responding oscillator lengths are l z and Zq with l z <C Zq. 
Throughout this Letter, we consider dipoles with no in- 
ternal spin degrees of freedom. 

The effective interaction V^i^r, 0) in the xy-plane is 
obtained by integrating out the harmonic motion in the 
z-direction. This yields 
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where £ = r 2 / (2Z 2 ) and Kq and K\ are irregular modified 
Bessel functions. (For the special cases = 90° and 9 = 
0°, expressions are also given in [12] ITS].) For electric 
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magnetic dipoles. The effective interaction ([2| reduces 
to the bare interaction dl]) when r ^> l z . Depending 
on the angle 9, V2D (r, (p) contains both attractive and 
repulsive regions. For 9 = 90° it is isotropic in the xy- 
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plane, whereas it is anisotropic with increasing attractive 
regions for decreasing angle. Since the focus is on Wigner 
states in this Letter, we restrict our analysis to the case 
arccos A= ~ 54.7° < < 90°, for which the interaction 
without the ^-integration is purely repulsive. 

Our numerical results are obtained by exact diagonal- 
ization using a basis of 2D harmonic oscillator orbitals. 
The multiparticle basis space is truncated by including 
all states with (kinetic and potential) energy up to some 
cutoff value E cut . Here, good convergence is achieved at 
E cut < 20Huj : and variations of E cut do not qualitatively 
affect our conclusions. Because of the computational ef- 
fort involved, the method is limited to particle numbers 
TV < 4. In the following we set l z /lo = 0.1. Our re- 
sults are not sensitive to this ratio as long as it is suf- 
ficiently small. In Fig. [2] we plot the particle density of 
the ground state for two different tilt angles O for a sys- 
tem with N = 3 dipolar bosons and fermions. For weak 
coupling, D = 0.1, the dipoles essentially form an ideal 
gas. The dip in the density for the fermions in Fig. [2] 
is simply a shell effect. With increasing interaction, the 
dipoles localize in Wigner states. Comparing the inter- 
action energy with the confinement kinetic energy yields 
the condition D 2 ^> 1 for Wigner crystallization. Our 
numerical results confirm this showing that the crystal- 
lized structure emerges continuously for D > 1 both for 
bosons and fermions for this finite-size system; for D = 5 
the localization can be clearly seen in Fig. [2] 

The Hamiltonian and the density of the ground state 
have rotational symmetry for O = 90°. To demon- 
strate particle localization, we therefore study the pair- 
correlated density p(r,r') = (r)^ ^(r / )'0(r / )^(r)) . 
Fixing one particle at position r, it gives the probabil- 
ity distribution of the remaining N — 1 particles. We see 
from Fig. [2] that for strong interactions and O = 90°, 
the particles localize at the vertices of an equilateral tri- 
angle. For O = 55°, the interaction is anisotropic in 
the xy-plane and almost vanishing between two dipoles 
whose relative vector is parallel to the x-axis. Figure [2] 
shows that the ground state for D = 5 features the three 
particles on a straight line for both bosons and fermions. 
This exploits the weaker regions of the interaction and 
minimizes the energy. A classical analysis minimizing 
the trap and interaction energies of three dipoles in a 2D 
trap indeed yields that the lowest energy configuration is 
all particles on a line along the x-axis for O < 62° and 
the three particles in an equilateral triangle for O > 62°. 

We now examine the transition from a line shaped to 
a triangular shaped ground state with increasing angle. 
Figure [3] depicts the lowest energy states of the three 
dipoles and their density as a function of the tilt angle 
O. In agreement with the classical analysis, the ground 
state changes from a linear to a triangular Wigner state 
as O is increased. For D = 5, the transition occurs at 
O ~ 56° for fermions and O ~ 57° for bosons which is 
somewhat lower than the classical prediction O ~ 62° due 
to quantum fluctuations. The density of the triangular 
state appears rotationally symmetric since the energy of 




Figure 2: Particle densities for three different coupling 
strengths, with dipole tilt angles G = 90° and G = 55°. The 
insets show the pair-correlated density functions; the crosses 
mark the position of one particle. The x and y axes range 
from — 4/ to 4/ , both here and in Figs. [3] and [4] 



three aligned classical dipoles in an equilateral triangle 
can be shown to be independent of its orientation. 

There is an intriguing difference between the transi- 
tions from linear- to triangular-shaped Wigner states for 
fermionic and bosonic dipoles: For fermions, the transi- 
tion is a sharp crossing, whereas, for bosons, it is a con- 
tinuous mixing of the triangular and linear states corre- 
sponding to an avoided crossing. This can be understood 
by considering that the Hamiltonian is invariant under 
the mirror operations P x : (x,y) —> (—x,y) and P y : 
(x,y) —> (x,—y). The energy eigenstates must there- 
fore also be eigenstates of P x and P y with eigenvalues 
±1. The key point is that for the Wigner states, one 
can infer their eigenvalues under these symmetry opera- 
tions from their quantum statistics. Consider the three- 
particle line configuration discussed above. Applying P x 
simply corresponds to exchanging the coordinates of the 
two outer localized particles. For bosons, this must give 
a plus sign due to the symmetrization of the wave func- 
tion. For fermions one must correspondingly get a mi- 
nus sign. Likewise, since all three particles have a finite 
probability to appear on the x-axis, the eigenvalue for P y 
must be +1 for both bosons and fermions. We write this 
as (P x , P y )B,iine = (+? +) f° r the bosonic line state and 
(P x , P 2/ )F,iine = ( — > +) f° r the fermionic line state. 

The triangular Wigner state is more subtle, as it con- 
tains components with the triangle in all possible orien- 
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Figure 3: The four lowest energy states for three strongly 
interacting bosonic (top panel) and fermionic (bottom panel) 
dipoles. By changing the dipole tilt angle B, the ground state 
changes from a Wigner state with line geometry to one with 
triangular geometry. The signs in brackets denote x- and y- 
parities, as explained in the text. 



tat ions. Take the component with one of the vertices 
of the triangle lying on the x axis: Operating with P y 
on this component corresponds to swapping the particles 
located at the two other vertices of the triangle symmet- 
rically placed above and below the x-axis; for bosons and 
fermions this will give the eigenvalue +1 and —1 respec- 
tively for this component. Since the Wigner state has 
to be an eigenfunction of P y , all other components of 
the wave function corresponding to tilted triangles must 
come in symmetric/antisymmetric pairs mirrored in the 
x-axis, so that the total wavefunction has the eigenvalue 
+1 (bosons) or —1 (fermions) for P y . The same anal- 
ysis applies to P x and we conclude that for the trian- 
gular Wigner state (P x , P y ) B , triangle = (+, +) for bosons 
and (P a ,,P 2/ )F,trian g ie = (-, -) for fermions. Since energy 
levels corresponding to wave functions of the same sym- 
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Figure 4: (a) The four lowest energy states for very strong 
interaction (axes as in Fig. 3). (b) Density contours for N = 4 
dipolar fermions. 



metries cannot cross, the arguments above explain why 
the line-triangle transition for bosons corresponds to an 
avoided crossing, while it is sharp for fermions. It is per- 
haps surprising that quantum statistics plays an impor- 
tant role for the transitions between essentially classical 
Wigner states. This can, however, be understood by not- 
ing that the particles can overlap during the transitions 
meaning that quantum statistics matters. 

It is illuminating to connect the states in Fig. [3] to the 
eigenstates of the angular momentum L z along the z-axis 
for = 90°: The lowest level in Fig. [3] naturally develop 
into to L z = for = 90°, whereas the next two levels 
correspond to L z ± 3 and not L z = ±1 as one would per- 
haps expect. This is because the strong repulsion favors 
states with large \m\ which are more spatially extended. 

In Fig. [4] (a), we plot the four lowest energy states in 
the limit of very strong interaction D = 10 (the interac- 
tion strength scales as D 2 ). Comparing with Fig. [3J we 
see that the difference between the bosons and fermions 
decreases with increasing coupling. This is as expected 
since the system approaches classical behavior in this 
limit. Reflecting this, the transition regions between dif- 
ferent ground states where quantum statistics is impor- 
tant are smaller for D = 10 than for D = 5: The size of 
the discussed anticrossing is AE = 0.53 for D = 5 and 
AE = 0.22 for D = 10. Also, the transition region has 
moved to larger angle c± 59° approaching the classical 
value ~ 62° with increasing D. 

In Fig. [4] (b) , we plot the density of the ground state 
of N = 4 fermionic dipoles for D = 5. In this case, 
the dipoles localize in a rhombic geometry that depends 
on the tilt angle. To examine the case of higher parti- 
cle numbers, we present in Fig. [5] the result of a classi- 
cal minimization of the interaction and trap energy for 
N = 19 dipoles. The ground state configurations were 
found using an iterative self-consistent procedure adjust- 
ing the particle positions to find the minimal energy. We 
see that in analogy with the N — 3, 4 cases considered 
above, the geometry of the lattice depends on the tilt 
angle. For = 90°, the particles localize in a hexago- 
nal lattice [13] . These classical calculations indicate that 
our main conclusions for quantum-mechanical few-body 
systems should be relevant also for larger systems. 

The strong coupling effects considered here are reach- 
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Figure 5: A classical calculation of the minimum energy con- 
figurations of N = 19 dipoles in a 2D trap. 



able using typical experimental numbers. The KRb 
molecule trapped by the JILA group [7] has an electric 
dipole moment of 0.57Debye whereas the RbCs molecule 
studied by the Innsbruck group has an electric dipole 
moment of 1.25Debye jT7} [18]. Trapping lengths of 
lo = 1 jam or Iq = 0.1 jam yield coupling strengths be- 
tween 0.6 < D 2 < 50. Also, there are alkali dimers with 



even larger dipole moments |6J. Using optical lattices, 
few-body systems with ultracold molecules may be real- 
ized, see e.g. the experiment in Ref. |4J. 

To summarize, we studied systems of strongly interact- 
ing bosonic and fermionic dipoles in a 2D harmonic trap, 
using both exact diagonalization and classical analysis. 
The dipolar interaction was shown to lead to a rich va- 
riety of Wigner states with nontrivial geometries which 
depend on the tilt angle of the dipoles with respect to the 
plane. Even though the Wigner states themselves can 
be well understood from a classical analysis, the tran- 
sitions between different geometries as the tilt angle is 
changed depends crucially on the quantum statistics of 
the dipoles. We showed how the effects analyzed here are 
relevant for typical experimental parameters. 
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